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ABSTRACT 

, Context. To understand the formation and evolution of galaxies, it is important to have a full comprehension of the role played by the 

^vq ■ metallicity, star formation rate (SFR), morphology, and color. The interplay of these parameters at different redshifts will substantially 

affect the evolution of galaxies and, as a consequence, the evolution of them will provide important clues and constraints on the galaxy 
evolution models. In this work we focus on the evolution of the SFR, metallicity of the gas, and morphology of galaxies at low redshift 
in search of signs of evolution. 

Aims. To analyze the S2N2 (logCHa/fS n]) vs. log(Ha/[Nn])) diagram as a possible segregator of star-forming, composite, and AGN 
galaxies, to study the evolution of the Baldwin, Phillips & Terlevich (1981) diagrams, as well as the evolution of the SFR, metallicity, 
and morphology, through the mass-metallicity, luminosity-metallicity, SFR-stellar mass, and SFR-metallicity relationships of star- 
forming galaxies from SDSS-DR5 (Sloan Digital Sky Survey-Data Release 5), using redshift intervals in bins of 0.1 from ~0 to 0.4. 
Methods. We used data processed with the STARLIGHT spectral synthesis code, correcting the fluxes for dust extinction, and esti- 
mating metallicities using the i? 2 3 method. We use the S2N2 diagnostic diagram as a tool to classify star-forming, composite, and 
AGN galaxies. We analyzed the evolution of the three principal BPT diagrams, estimating the SFR and specific SFR (SSFR) for our 
samples of galaxies, studying the luminosity and mass-metallicity relations, and analyzing the morphology of our sample of galaxies 
, through the g — r color, concentration index, and SSFR. 

Results. We found that the S2N2 is a reliable diagram to classify star-forming, composite, and AGNs galaxies. We demonstrate that 
the three principal BPT diagrams show an evolution toward higher values of [O in] /15007/H/3 due to a metallicity decrement. We 
found an evolution in the mass-metallicity relation of ~ 0.2 dex for the redshift range 0.3 < z < 0.4 compared to our local one. 
From the analysis of the evolution of the SFR and SSFR as a function of the stellar mass and metallicity, we discovered a group of 
galaxies with higher SFR and SSFR at all redshift samples, whose morphology is consistent with those of late-type galaxies. Finally, 
the comparison of our local (0.04 < z < 0.1) with our higher redshift sample (0.3 < z < 0.4), show that the metallicity, the SFR and 
morphology, evolve toward lower values of metallicity, higher SFRs, and late-type morphologies for the redshift range 0.3 < z < 0.4. 
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1 . Introduction burst line on the three BPT diagrams. Kauffmann et al. (2003a, 

hereafter Kauf03) shifted this starburst limit to a lower and 

Baldwin, Phillips & Terlevich (1981, hereafter BPT) were the more precise one in the [Nn]/Ha diagram, excluding Seyfert- 

nrst to propose diagnostic diagrams to classify galaxies into Hii composite objects whose spectra contain significant contri- 

$J ; starburst or active galactic nucleus (AGN), based on the dom- butions from both AGN and star formation, from pure SF galax- 

<P_ \ inant energy source in emission-line galaxies, since AGNs have i es . Galaxies between the Kauf03 and KewOl divisions are con- 

a much harder ionizing spectrum than hot stars. Revised and re- sidered as composite galaxies. 

fined by Veilleux & Osterbrock (1987), the three BPT empirical There ws more division criteria between SF galaxies and 
diagnostic diagrams use the optical line ratios [Oi] ^6300/Ha, AGNs in the [Nn]/H« vs [Om] A5007/H/3 BPT diagram, such 
[S n] AA6717, 6731/Ha, [Nn] ^6583/Ha, and [Om] A5001/R/3 as the ne of Stasiriska et al. (2006), which used a limit lower 
(hereafter [Nn] will refer to [Nn] ,16583, and [Sn] to [Sn] man that of K auf03, based on a more rigurous criterion, and 
AA6717, 6731). The BPT diagrams are the most widely used the one of Lee et al. (2007), which used an intermediate em- 
method to segregate between star-forming galaxies and AGNs, pirical line between the Kauf03 and KewOl divisions. It is pos- 
since the lines in star-forming (SF) galaxies are emitted by H n sib i e , however, to classify SF galaxies and AGNs using only the 
regions, which are ionized by massive stars, while AGNs are [Nn]/Ha ratio, as discussed in Stasiriska et al. (2006), since the 
ionized by a harder radiation field. l e f t arm of the [N n ]/H a diagram (see for example Fig. 2) is a 
Kewleyetal. (2001, hereafter KewOl) used a combination of measure of the combination of the metallicity and the ioniza- 
stellar population synthesis models and detailed self-consistent tion parameter. Then, larger values of this ratio indicate that the 
photoionization models to create a theoretical maximum star- galaxy host an AGN. Stasiriska et al. (2006) classify as SF galax- 
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ies those with log([Nn]/Ha) < -0.4, composite galaxies those 
with -0.4 < log([Nn]/Har) < -0.2, and AGNs those galaxies with 
log([Nn]/Hor) > -0.2. 

Following with the objective of segregate SF from compos- 
ite and AGNs galaxies, in this work we study the S2N2 dia- 
gram as a reliable segregator of galaxies. This log(Ha/[S n]) vs. 
log(Ha/[Nn]) diagram was introduced by Sabbadin et al. (1977) 
as a useful tool to separate galactic planetary nebula (PNe), 
H ii regions, and supernova remnants (SNRs). This diagram was 
later applied to Herbig-Haro objects (Canto 1981), Galactic PNe 
(Garcfa-Lario et al. 1991, Riesgo & Lopez 2005), and extra- 
galactic PNe (Magrini et al. 2003). The S2N2 diagram has been 
used also as a metallicity and ionization parameter indicator for 
extragalactic H n regions by Viironen et al. (2007). 

The S2N2 diagram has been also applied to galaxies by some 
authors. For example, Moustakas & Kennicutt (2006) studied 
whether there was a difference between integrated spectra of 
galaxies and the spectra of individual Hn regions. Dopita et al. 
(2006) used the S2N2 diagram, among others, for abundance di- 
agnostics using photoionization models. Nevertheless, the [S n] 
flux shows always deficiences when generated by photoioniza- 
tion models (e.g. Levesque et al. 2010). Also, Lamareille et al. 
(2009) and Perez-Montero et al. (2009) used the S2N2 diagram 
as a segregator of SF from Seyfert 2 galaxies, but using different 
ratios: log([Nn]/Ha) vs log([S n]/Ha). However, in their divi- 
sion Lamareille et al. (2009) do not distinguish between SF and 
composite galaxies, also, they used equivalent widths instead of 
emission line fluxes, which could affect the results (Kobulnicky 
& Kewley 2004). 

The formation and evolution of galaxies at different cosmo- 
logical epochs are driven mainly by two linked processes: the 
star formation history and the metal enrichment. Thus, from an 
observational point of view, the star formation rate (SFR), the 
metallicity and the stellar mass of the galaxies at different epochs 
will give us important clues on the evolution of galaxies. The 
first quantitative SFRs were derived from evolutionary synthe- 
sis models of galaxy colors (Tinsley 1968, 1972, Searle et al. 
1973), confirming the trends in SFRs and star formation histo- 
ries along the Hubble sequence, and giving the first predictions 
of the evolution of the SFR with cosmic lookback time. The de- 
velopment of more precise direct SFR diagnostics includes the 
integrated emission-line fluxes (Cohen 1976, Kennicutt 1983), 
near-ultraviolet continuum fluxes (Donas & Deharveng 1984), 
and infrared continuum fluxes (Harper & Low 1973, Rieke & 
Lebofsky 1978, Telesco & Harper 1980); see Kennicutt (1998) 
for a review. The hydrogen Balmer line Ha is currently the most 
reliable tracer of star formation, since in H n regions and star- 
forming galaxies, the Balmer emission-line luminosity scales 
directly with the total ionizing flux of the embedded stars. A 
widely known calibration of the Ha line as SFR tracer is the 
one devised by Kennicutt (1998). However, it is important to 
take into account corrections for stellar absorption and redden- 
ing to obtain SFRs in agreement with the ones derived using 
other wavelengths (e.g. Rosa-Gonzalez et al. 2002, Chariot et al. 
2002, Dopita et al. 2002). In parallel, other diagnostics have been 
developed using the oxygen doublet [O n] /13726, 3729 for the 
redshift range z ~ 0.4- 1.5 (e.g. Gallagher et al. 1989, Kennicutt 
1998, Rosa-Gonzalez et al. 2002, Kewley et al. 2004). Moreover, 
this diagnostic is usefull when the Ha line is not easily observ- 
able at higher redshifts (z > 0.4 in the optical). However, the 
[O n] doublet presents problems in reddening and abundance de- 
pendence (Jansen, Franx & Fabricant 2001, Chariot et al. 2002). 
Alternatively, it is possible to estimate the SFR from the soft X- 



ray luminosity, which is comparable to that determined from the 
Ha luminosity (Rosa Gonzalez et al. 2009, Rovilos et al. 2009). 

A strong dependence of the SFR and the stellar mass and its 
evolution with redshift has been found, with the bulk of star for- 
mation occurring first in massive galaxies, and later in less mas- 
sive systems (e.g. Guzman et al. 1997, Brinchmann & Ellis 2000, 
Juneau et al. 2005, Bauer et al. 2005, Bell et al. 2005, Perez- 
Gonzalez et al. 2005, Feulner et al. 2005, Papovich et al. 2006, 
Caputi et al. 2006, Reddy et al 2006, Erb et al. 2006, Noeske et 
al. 2007a, Buat et al. 2008). In the local universe, several stud- 
ies have illustrated a relationship between the SFR and stellar 
mass, identifying two populations: galaxies on a star-forming se- 
quence, and "quenched" galaxies, with little or no detectable star 
formation (Brinchmann et al. 2004, Salim et al. 2005, Lee 2006). 
At higher redshift, Noeske et al. (2007a) showed the existence of 
a "main sequence" (MS) for SF galaxies in the SFR-stellar mass 
relation over the redshift range 0.2 < z < 1.1. From the galaxies 
considered in this study, it was shown that the slope of the MS 
remains constant to z > 1, while the MS as a whole moves to 
higher SFR as z increases. 

Metallicity is another important property of galaxies, and its 
study is crucial for a deep understanding of galaxy formation 
and evolution, since it is related to the whole past history of the 
galaxy. Metallicity is a tracer of the fraction of baryonic mass 
already converted into stars and is sensitive to the metal losses 
due to stellar winds, supernovae and active nuclei feedbacks. A 
detailed description of the different metallicity methods and cal- 
ibrations are given in Lara-Lopez et al. (2009a,b). 

Stellar mass and metallicity are strongly correlated in SF 
galaxies, with massive galaxies showing higher metallicities 
than less massive galaxies. This relationship provides crucial in- 
sight into galaxy formation and evolution. The mass-metallicity 
(M-Z) relation was first observed by Lequeuxet al. (1979), has 
been intensively studied (Skillman et al. 1989; Brodie & Huchra 
1991; Zaritsky et al. 1994; Richer & McCall 1995; Garnett et 
al. 1997; Pilyugin & Ferrini 2000, among others), and it is well 
established by the work of Tremonti et al. (2004, hereafter T04) 
for the local universe (z ~ 0.1) using SDSS data. The study of 
the redshift evolution of the M - Z relation has provided us with 
crucial information on the cosmic evolution of star formation. 

Regarding the evolution of the M-Z relation for SF galax- 
ies at z < 1, Savaglio et al. (2005), have investigated the mass- 
metallicity relations using galaxies at 0.4 < z < 1, finding that 
metallicity is lower at higher redshift by ~ 0.15 dex. Moreover, 
Maier at al. (2005), Hammer et al. (2005), and Liang et al. (2006) 
found that emission line galaxies were poorer in metals at z ~ 
0.7 than present-day spirals. A study of Lamareille et al. (2009) 
focused on the evolution of the M - Z relation up to z ~ 0.9, 
suggesting that the M-Z relation is flatter at higher redshifts. 
However, Carollo & Lilly (2001), from emission-line ratios of 
15 galaxies in a range of 0.5 < z < 1, found that their metallic- 
ities appear to be remarkably similar to those of local galaxies 
selected with the same criteria. Also, Lilly et al. (2003), from a 
sample of 66 SF galaxies with 0.47 < z < 0.92, claim a smaller 
variation in metallicity of ~ 0.08 dex compared with the metal- 
licity observed locally, showing only modest evolutionary effects 
(for more details about the M - Z relation, see Lara-Lopez et al. 
2009b). 

In a recent study, Calura et al. (2009) have demonstrated the 
importance on the morphology of galaxies when deriving the 
M-Z relation since, at any redshift, elliptical galaxies present 
the highest stellar masses and the highest metallicities, whereas 
the irregulars are the least massive galaxies, characterised by the 
lowest O abundances. 
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In this paper, we consistently approach several topics, start- 
ing with the introduction of the S2N2 as a reliable diagram to 
classify galaxies, the analysis of the metallicity evolution of 
galaxies in the three BPT diagrams, and, for a better understand- 
ing of the processes involved in the observed evolution of galax- 
ies at low redshift, we studied the mass, metallicity and SFR 
relations, such as the M-Z, metallicity-SFR and mass-SFR rela- 
tions. We also point out that the morphology of galaxies play an 
important role when deriving conclusions, since late-type galax- 
ies will result in lower metallicity estimates and higher SFRs 
than early-type (Calura et al. 2009). 

This paper is structured as follows, in Sect. 2 we detail the 
data used for this study, the dust extinction correction and the 
metallicity estimates for our sample of galaxies, in Sect. 3 we 
introduce the S2N2 as a reliable diagram to segregate SF, com- 
posite, and AGNs galaxies. In Sect. 4 we analyzed the evolu- 
tion of the BPT diagrams. In Sect. 5 we investigate the evolu- 
tion of the mass-metallicity and luminosity-metallicity relations, 
whereas in Sect. 6 we discuss the relations between the SFR 
and SSFR with stellar mass and metallicity, as well as the mor- 
phology of our galaxies using colors, concentration index, and 
SSFRs. Finally, conclusions are given in Sect. 7. 

2. Data processing and sample selection 

We selected emission line galaxies from SDSS-DR5 (Adelman- 
McCarthy et al. 2007). Data were taken with a 2.5 m telescope 
located at Apache Point Observatory (Gunn et al. 2006). The 
SDSS spectra were obtained using 3 arcsec diameter fibres, cov- 
ering a wavelength range of 3800-9200 A, and with a mean 
spectral resolution A/AA ~ 1800. The SDSS-DR5 spectroscopy 
database contains spectra for ~ 10 6 objects over ~ 5700 deg . 
Further technical details can be found in Stoughton et al. (2002). 

We used the SDSS-DR5 spectra from the STARLIGHT 
databaseQ, which were processed with the STARLIGHT spectral 
synthesis code, developed by Cid Fernandes and collaborators 
(Cid Fernandes et al. 2005, 2007, Mateus et al. 2006, Asari et al. 
2007). From the spectra, the STARLIGHT code subtracts the 
continuum, obtaining the emission lines fluxes measurements 
for each galaxy. For each emission line, the STARLIGHT code 
returns the rest frame flux and its associated equivalent width, 
linewidth, velocity displacement relative to the rest-frame wave- 
length and the S/N of the fit. In the case of Balmer lines, the un- 
derlying stellar absorption was corrected by the STARLIGHT 
code using synthetic spectra obtained by fitting an observed 
spectrum with a combination of 150 simple stellar populations 
(SSPs) from the evolutionary synthesis models of Bruzual & 
Chariot (2003). 

From the full set of galaxies, we only consider galaxies 
whose spectra show in emission the Ha, H/3, [Nn], [On] A3121, 
[Om] /14959, [Om] /15007, [Oi] ,16300, [Sn] lines. We selected 
galaxies with a signal-to-noise ratio higher than 3cr for the Ha, 
Hy6, and [N n] lines. 

In order to identify any evolution of galaxy parameters or 
relations, we divided our sample in four redshift intervals as 
follow: 0.04 < zo < 0.1, 0.1 < z\ < 0.2, 0.2 < z 2 < 0.3, 
0.3 < Z3 < 0.4. The lower limit of zo corresponds to an aper- 
ture covering fraction of 20%, which is the minimum required 
to avoid domination of the spectrum by aperture effects (Kewley 
et al. 2005). This classification give us 85931 galaxies for zo, 
48888 galaxies forzi, 3278 galaxies forz2, and 199 galaxies for 

Z3- 

1 http://www.starlight.ufsc.br 



We selected galaxies with an apparent Petrosian r magni- 
tude of 14.5 < r < \1 .11 in the redshift samples zo, Z\, and 
Z2, which yields 82884, 44763, and 1802 galaxies, respectively, 
corresponding to the magnitude completeness at these redshifts 
(see Fig. 1). Galaxies of the Z3 sample have a different complete- 
ness range 16.9 < r < 18.8, as observed in Fig. 1, giving 119 
galaxies. We used the zo and z\ sample of galaxies with its re- 
spective completeness, but for galaxies of samples zi and zi we 
used both, those in the completeness range, and those out of the 
completeness range. The reason for this, is to improve the galaxy 
statistics by increasing the number of them. As we shown in the 
next sections, the main results are similar using galaxies in the 
magnitude completeness and galaxies of the total sample. 




14 16 18 20 

Fig. 1. Normalized histogram of the apparent Petrosian r magni- 
tudes in the four redshift bins. Dark solid line represents galaxies 
at zo, dashed line the galaxies at zi, dot-dashed lines the galax- 
ies at Z2, and clear solid line, the galaxies at Z3- The black arrow 
shows the completeness limit for the samples zo, Z\, and zz, and 
clear arrow points the same for zi- 



2. 1. Sample selection for Section 3 

In Section 3, to study the S2N2 diagram as a segregatorof differ- 
ent types of galaxies, we adopt the "main galaxy sample" (e.g. 
Strauss et al. 2002) with Petrosian r magnitudes in the range 
14.5 < r < 17 .77, and the redshift interval zo, taking into ac- 
count all the emission lines and the signal-to-noise ratio men- 
tioned above, which yields 82884 galaxies. In this Section, we 
used AGN, composite, and SF galaxies. 

2.2. Sample selection for Section 4 

In Section 4, we study the evolution of the three BPT diagrams 
using the four redshift intervals, magnitude intervals, and the 
signal-to-noise restrictions mentioned above. In this Section, we 
used AGN, composite, and SF galaxies. Also, to study the metal- 
licity evolution of the SF galaxies of the BPT diagrams (see Fig. 
6), we used the sample of Sections 5 and 6, mentioned below. 
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2.3. Sample selection for Sections 5 and 6 

For Sections 5 and 6, we selected SF galaxies following the 
criterion given by Kauf03 in the BPT empirical diagnostic di- 
agram: log[Om] /15OO7/H V < 0.61/{log([Nn] /Ha)-0.05} + 1.3, 
the same used by Veilleux & Osterbrock (1987), Kewley et al. 
(2001, 2006), and Stasiriska et al. (2006), among others. After 
all these selections, the number of galaxies of each redshift bin 
is reduced to 61921 SF galaxies for zo, 27853 for z\ , 1671 for zi, 
and 67 H n galaxies for Z3. 

The extinction correction and metallicity estimates were cal- 
culated as in Lara-Lopez et al. (2009b). Our extinction correc- 
tion was derived using the Balmer decrements in order to ob- 
tain the reddening coefficient C(H/T). We used the Cardelli et al 
(1989) law, with R v = A,,/E(B - V)=3.1, assuming case B recom- 
bination with a density of 100 crrT 3 and a temperature of 10 4 K, 
with HafH/3 = 2.86 (Osterbrock 1989). 

We estimated metallicities using the R23 relation intro- 
duced by Pagel et al. (1979), fl 2 3=([On] ^3727+[Ora] AA4959, 
5007)/H/3, and adopted the calibration given by Tremonti et al. 
(2004), 12 + log(0/H) = 9.185 - 0.313x - 0.264x 2 - 0.321x 3 , 
where x = log 7?23- We selected the upper branch of the double- 
valued R23, in which the Tremonti et al. (2004) calibration is 
valid, taking 12+log(0/H) > 8.4 and log([Nn]/[On]i3727) > - 
1.2, since the upper and lower branches of the R23 calibration 
bifurcates at those values (see Kewley & Ellison 2008). 

Applying this final criterion, we end with 58866 galaxies for 
Zo, 24385 forzi, 1631 forz2 (from which 712 galaxies are in their 
completeness magnitude interval), and 62 galaxies for Z3 (from 
which 41 galaxies are in their completeness magnitude interval), 
all of them in the upper branch of the R23 relation, corresponding 
to the ~99 % to the Hn classified galaxies. Then, we are not 
introducing a bias selecting the upper branch for this samples. 



3. The S2N2 diagnostic diagram as a star-forming, 
composite, and AGN galaxies segregator 

As mentioned in Section 1, BPT diagrams are the most used 
method to segregate between star-forming and AGN galaxies. 
From the three BPT diagrams ([Nnl/Ha, [Sn]/Ho', and [Oi] 
^6300/Hor vs [O m] ^5007/HyS), the most used is the [N n]/Ua vs 
[O m] /15007/H/3 one, since it is the only one that can segregate 
pure SF and composite galaxies, as demonstrated by Kewley et 
al. (2006) and Perez-Montero et al. (2009). The other two BPT 
diagrams are not useful for segregating SF from composite ob- 
jects. 

Commonly known as the S2N2, the log(Ha/[S 11]) vs. 
log(Ho'/[Nii]) diagram has been used to separate planetary neb- 
ulae (PNe), Hn regions, and supernova remnants (SNRs, see 
Sabbadin et al 1977, Riesgo & Lopez 2006, Viironen et al. 
2007). We propose the S2N2 diagram to classify SF, compos- 
ite, and AGN galaxies, something until now only possible with 
the [Nn]/Ho' vs [Om] A5007/H/3 diagram. However, the S2N2 
diagram use only the Ha, [Nn], and [S 11] emission lines, all of 
them close in wavelength, avoiding reddening corrections, and 
making possible its use for surveys limited in spectral range. 

From the "main galaxy sample", we consider SF galaxies as 
those lying below the Kauf03 division, composite galaxies as 
those lying between the Kauf03 and KewOl lines, and AGNs 
galaxies, those above the KewOl division (see Fig. 2). From the 
total sample (82884 galaxies, see Sect. 2), the 71.4%, 19%, and 
9.6% correspond to SF, composite, and AGN galaxies, respec- 
tively. Taking this classification as a reference, we plotted the 
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Fig. 2. Log([Nn]/Ha) vs log([Om] A5007/H/3) BPT diagram. 
Solid line show the Kauf03 empirical division between SF and 
composite galaxies, and dashed line represents the KewOl star- 
burst limit. [See the electronic edition of the Journal for a color 
version of this figure.] 



S2N2 diagram with the three classifications of galaxies (see Fig 
3 and 4). 




-0.5 0.0 0.5 1.0-0.5 0.0 0.5 1.0 

log(Ha/[SII]U6717,6731) log(Ha/[SII]?A671 7,6731) 

Fig. 3. Contour plots for star-forming, composite, and AGN 
galaxies, (a) Star-forming, darker (blue) contour, and compos- 
ite, lighter (orange) contour, enclosing the ~90% (~1.8 cr) of 
each sample. The solid line establish our limit for star-forming 
galaxies, while dot-dashed line delimit the almost pure star- 
forming galaxies, (b) Composite, darker (orange) contour, and 
AGN, lighter (grey) contour, enclosing ~75% (~1.2 cr) of each 
sample. The dashed line shows our separation for composite and 
AGN galaxies, solid and dot-dashed lines as in panel (a). [See 
the electronic edition of the Journal for a color version of this 
figure.] 



In order to establish division lines to separate SF, composite 
and AGN galaxies in the S2N2 diagram, we generated contour 
plots for each category of galaxies (see Fig. 3). As our sample 
of galaxies is larger for the SF and composite galaxies, we 
used contours enclosing ~90% for those galaxies. However, as 
AGN galaxies are less numerous, we used contours enclosing 
~75% for composite and AGNs. The contour plots shown in 
Fig. 3a delimit two tangent parallel lines, generating with this 
criterium division lines defined by Eqs. (1) and (2). In order to 
define a division line between composite and AGN galaxies, we 
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Fig. 4. (a) log(Ha'/[S n] ) vs. log(Ha/[Nii]) diagram, corresponding to the star-forming galaxies below the Kauf03 separation in 
the BPT diagram. We overploted the grid of the photoionization models of Kewley et al. (2001) created with PEGASE, taking an 
instantaneous burst model with an electronic density of 10 cm -3 . Dot-dashed, solid, and dashed line, represent our division for pure 
star-forming and composite, star-forming and composite, and composite and AGN galaxies, respectively, (b) Same diagram with 
the grid of the photoionization models of Kewley et al. (2001) created with PEGASE, taking a continuous model with an electronic 
density of 10 cirT 3 . [See the electronic edition of the Journal for a color version of this figure.] 



sampled the plot area with parallel lines of Eq. (1) in bins of 
0.02 dex, generating in this way histograms for composite and 
AGN galaxies, where Eq. 3 corresponds to the intersection of 
both histograms. 

(i) Pure star-forming galaxies are separated by: 

log(Hor / [N n] ) > -0.28 x log(Har/ [S n]) + 0.5 1, ( 1 ) 

98.8% of galaxies above this line are SF galaxies and correspond 
to 88% of the SF sample. 

(ii) Star-forming and composite galaxies are divided by: 

log(Ha/[Nn]) > -0.28 x log(Har/[S n]) + 0.40, (2) 

99.8% of the SF sample lay above this line. However, from 
all galaxies above this line, 8% will correspond to composite 
galaxies, and the remaining 91% to SF galaxies. 

(iii) Composite and AGN galaxies are divided by: 

log(Ha/[Nn]) > -0.28 x log(Hor/[Sn]) + 0.25, (3) 

74.6% of the galaxies below this line are AGNs. Unfortunately, 
this diagram does not allow separating Seyfert from Liner galax- 
ies (see Fig. 3). Composite galaxies are selected as those be- 
tween the lines of Eq. (1) and (3). From all galaxies between 
both lines, 61.2%, 33.2%, and 5.6% correspond to composite, 
SF, and AGN galaxies, respectively. 

In table 1 we compare the N2 ratio, considering the Stasihska 
limits, with the S2N2 diagram in order to test their respective 
ability for segregating SF, composite, and AGN galaxies. This 
table shows the percentage of galaxies classified according to 
each division line, as well as its corresponding contamination 
and percentage of missed galaxies taking as a reference the clas- 
sification of the log([Nn]/Ha) vs log([Om] /15007/H/3) diagram. 



As can be appreciated for the SF segregation, the N2 division 
miss the highest percentage of SF galaxies, while our pure star- 
forming division (Eq. 1), although with the same contamination 
from composite galaxies, miss less SF galaxies. Moreover, our 
SF limit (Eq. 2) enclose almost all SF galaxies with a small 
contamination from composite galaxies. This provides the ad- 
ditional advantage with respect to the N2 diagram, that the user 
can choose between the possibilities of selecting either most SF 
galaxies, or galaxies with the smallest contamination from com- 
posite galaxies. The composite division for the N2 ratio shows 
a similar percentage of composite and SF galaxies. However, 
the S2N2 division, allows obtaining a higher percentage of com- 
posite galaxies, missing less composite galaxies, with a lower 
contamination of SF galaxies, and with only a quite small incre- 
ment of AGN contamination. For segregating AGNs, while the 
S2N2 diagram provides a lower contamination from composite 
galaxies, miss ~6% more AGNs than using the N2 ratio. 

In Fig. 4, we overplot the pre-run photoionization grids 
of Dopita et al. (2000) and KewOl for an instantaneous burst 
model and for a continuous starburst model. The best grid to 
our SF galaxies for the S2N2 diagram is the corresponding to 
an instantaneous burst model, with an electronic density of 10 
cm" 3 , using the PEGASE code (see Fig. 4a). As explained in 
Dopita et al. (2000), the high surface brightness isolated extra- 
galactic H ii regions are in general excited by young clusters of 
OB stars and that, in this case, the ionizing EUV spectra and 
H ii region emission-line spectra predicted by the PEGASE and 
STARBURST99 codes for an instantaneous, zero-age star for- 
mation model, are essentially identical. 

In their work, Dopita et al. (2000) and KewOl mod- 
eled a large sample of infrared starburst galaxies using both 
the PEGASE v2.0 (Fioc & Rocca-Volmerange 1977) and the 
STARBURST99 (Leitherer et al. 1999) codes to generate the 
spectral energy distribution (SED). In both cases, MAPPINGS 
III code was used to compute photoionization models. The pre- 
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34 


S2N2(1&3) 
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1.2 
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33.2 
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15.3 


25.4 


S2N2(2) 


91 
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8 

















Table 1. Comparison between the S2N2 diagram and the N2 ratio to segregate SF, composite, and AGN galaxies. S2N2(1&3) and 
S2N2(2) refer to galaxies taking as a reference Eq. (1&3) and Eq. (2), respectively. The SF block indicate for the S2N2(2) division, 
for example, that taking all galaxies above Eq. (2), the 91% will correspond to SF galaxies, and the rest 8% to composite galaxies, 
also, that we have missed 0.2% of the original SF sample. In the case of the composite block, for the S2N2(1&3) division we are 
taking galaxies between the lines defined by Eqs. (1) and (3). 



run grids use photoionization models with ionization parameters 
q (cm s~') in the range 5xl0 6 < q < 3xl0 8 , and metallicities 
from Z=0.05 to 3 Z Q ; moreover, two values for electronic den- 
sity were used, 10 and 350 crrT 3 . 

For starburst galaxies, it is expected to have a continuous star 
formation over at least a galactic dynamical timescale, then, the 
assumption of a continuous rather than an instantaneous burst 
of star formation would be more accurate. We generated the 
PEGASE grid of KewOl for a continuous starburst model and 
found, as expected, that this corresponds to SF and composite 
objects (see Fig. 4b), since this is the limit used by KewOl to 
parametrize an extreme starburst line in the BPT diagrams. 

The photoionization grids generated with the 
STARBURST99 code were not hard enough to produce 
the needed [Sn] flux to enclose all galaxies in the S2N2 
diagram. The BPT diagrams are most sensitive to the spectral 
index of the ionizing radiation field in the 1-4 ryd interval, and 
the PEGASE ionizing stellar continuum is harder in this range 
than that of STARBURST99, being PEGASE the only models 
that encompass nearly all of the observed starburst on all three 
of the BPT diagrams. We also tried the grids of Levesque et al. 
(2010), which use the STARBURST99 code, but these grids 
comprise an insufficiently hard ionizing radiation field, leading 
to deficiencies in the [S n] fluxes produced by the models. 

In galaxy surveys there are at least three methods commonly 
used to deal with the presence of AGNs. The first one consists on 
removing galaxies hosting AGNs by cross-correlating the sam- 
ple with published AGN catalogues (e.g. Condon et al. 2002, 
Serjeant et al. 2002). The second one deals with the identification 
of the galaxies through the so called BPT diagram. Nevertheless, 
at z > 0.5, the Ha line is redshifted out of the optical range. The 
third method consists in subtracting AGNs in a statistical man- 
ner, used when no other methods are applicable (e.g. Tresse & 
Maddox 1998). 

In addition to those methods, we propose the use of the S2N2 
diagram, which has demonstrated the ability of accurately seg- 
regate SF from composite and AGN galaxies. The S2N2 dia- 
gram has the following advantages: it is not necessary an extinc- 
tion correction, since all the emission lines are close in wave- 
length; requires only a small spectral range, making it suitable 
for surveys of limited spectral coverage; the SF and composite 
divisions of the S2N2 diagram offer less contamination in all 
cases, with respectively higher number of galaxies for SF and 
composites than using only the N2 ratio; additionally, the user 
can choose any of the SF galaxies divisions provided for the 
S2N2 diagram, either if the smallest contamination from com- 
posite galaxies (Eq. 1), or selecting the most SF galaxies (Eq. 2) 
is required. Finally, the use of the [S n] lines do not reduce the 
number of galaxies, since those lines for AGN galaxies are, in 
mean ~ 1 .3 times stronger than the [O m] line used in the BPT di- 



agrams. Then, comparing the number of galaxies of all the types, 
the S2N2 diagram have 1% more galaxies than the log([Nn] 
/Ha) vs log([0 m] /HjS) diagram. Although this diagram has been 
used in the past for galaxies, it is the first time that it is presented 
as a diagnostic diagram for classifying galaxies. Given its ad- 
vantages, we propose the use of this diagram as an alternative to 
the BPT diagrams and N2 ratio to classify galaxies. 



4. Evolutionary effects on the BPT diagrams 

As explained in Section 3, the BPT and other optical emission 
lines diagnostic diagrams have become important in the classi- 
fication of galaxies. In this section, our aim is to investigate the 
effects of the evolution of galaxies from the three BPT diagrams. 
For this purpose, and with the objective of increase our number 
of galaxies, we did not take any restriction in magnitude, as de- 
tailed in the sample selection. 

In Fig. 5 we show the three BPT diagrams for the four 
redshift samples. As redshift increases, we observe that [Om] 
i5007/Hj6 goes toward higher values. In order to explain this 
shift, in Fig. 6 we plotted the ratio [Oiii] /i5007/Hj8 versus 
redshift and metallicity only for SF galaxies selected with the 
Kauf03 criterion. The gap observed around z ~ 0.145 (see Fig. 
6a, c) is due to the H/3 line falling nearby the 5577 A sky line, 
because the residuals are significant and, as a consequence, mea- 
surements of Hj6 around this redshift were lost. As shown in Fig. 
6b, there is a clear tendency of the [O m] /15007/H/? ratio towards 
higher values with redshift, which is explained by examining the 
same ratio against 12+log(0/H). The ratio [Oiii] i5007/H/3 has 
demonstrated to correlate linearly with metallicity (see, for ex- 
ample, Liang et al. 2006). Then, a decrement in 12+log(0/H) 
will result in higher values of [Om] /15007/H/3 (see Fig. 6b). 
We observe a decrement of ~0.2 dex in [O m] /15007/H/?, and a 
decrement of ~0.1 dex in 12+log(0/H) for the Z3 redshift range 
with respect to the zq range. 

In previous papers (Lara-Lopez et al. 2009a, b) we reported 
a decrement in 12+log(0/H) of ~0.1 dex for the redshift range 
0.3 < z < 0.4 comparing galaxies in the same range of luminos- 
ity at different redshift intervals. Since the possible bias, such 
as luminosity, mass and aperture effects of those samples were 
carefully studied, we demonstrated there that this decrement in 
metallicity is due to an intrinsic evolution of the galaxies. 

Although our Z3 sample corresponds to luminous galaxies, 
if we compare galaxies with the same luminosity, taking as a 
reference our previous papers, the metallicity decrement will be 
again of ~0.1 dex, and as consequence, the effects on the BPT 
diagrams will be the same. Therefore, the evolution observed in 
the [O m] /15007/H/? lines ratio toward higher values in the three 
BPT diagrams, could be attributed to a metallicity evolution. 
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log([NII]\6583/H«) log([SII]671 7,673 1 /Ha) log([OI]/Ha) 



Fig. 5. Evolution of line ratios from the BPT diagrams. Black dots show galaxies in the redshift range zo, blue dots galaxies in 
the Zi range, red big triangles represent galaxies in the completeness of the zi range, while small red triangles galaxies out of 
the completeness, yellow big circles represents galaxies in the completeness of the zj range, and small circles galaxies out of the 
completeness. The black solid line shows the Kauf03 limit for SF galaxies, and the red solid line shows the KewOl limit for starburst 
galaxies in the three BPT diagrams. 




Fig. 6. Evolution with redshift of the log([Ora] A5007/H/3) ratio for star-forming galaxies (left), and 12+log(0/H) vs log([Om] 
ASOOl/Hfi) for the same galaxies (right). Symbols are the same as in Fig. 5. The gap observed at z ~ 0.145 is due to the 5577 A sky 
line falling in the H/3 line, missing galaxies around this redshift. 



On the other hand, we analyze the [N ii]/Hq' ratio against red- widely studied since it is not severely affected by dust extinc- 
shift and metallicity (see Fig. 6c, d). The [Nn]/Hor ratio is also tion (see Pettini & Pagel 2004). Among the calibrations of the 
a metallicity index, commonly known as N2, and it has been N2 index, we have for example those of Raimann et al. (2000), 
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Denicolo et al. (2002), and Liang et al. (2006). In Fig. 6d we 
observe a clear increasing trend of metallicity following the in- 
crease of the N2 index up to 12+log(0/H) ~9.0. The galaxies 
with 12+log(0/H) >9 show a flattening and a slightly decrease 
of the N2 index with metallicity (see Fig. 6d). This trend was 
explained by Kewley et al. (2002) using photoionization mod- 
els as follows: when the secondary production of nitrogen dom- 
inates, at somewhat higher metallicity, the [Niil/Ha line ratio 
continues to increase, despite the decreasing electron tempera- 
ture. Eventually, at still higher metallicities, nitrogen becomes 
the dominant coolant in the nebula, and the electron temperature 
falls sufficiently to ensure that the nitrogen line weakens with 
increasing metallicity. Liang et al. (2006), using SDSS galaxies 
with redshifts 0.04 < z < 0.25, observed a small decrement of 
the N2 index against metallicity; this turnover of the N2 index is 
more evident for the higher redshifts zi and zj, in our sample (see 
Fig. 6d). The turnover of the N2 ratio will produce in the [Om] 
/15007/H/3 vs. [N ii]/Hck BPT diagram, the turnover zone around 
N2 0.4, which is more evident in a density plot (see Fig. 2). 

Regarding the two left BPT diagrams of Fig. 5b and c, since 
they share the ratio [Om] A5001/H/3, the evolutionary effects 
due to a decrement in metallicity will be the same as discussed 
above. The ratio [S n]/Ha has never been used before as a metal- 
licity indicator because it is far more sensitive to ionization than 
to metallicity (Liang et al. 2006). Moreover, it is double-valued 
with metallicity (see Fig. 4), whereas the ratio [OI]/Ha is not a 
metallicity indicator. 

Therefore, after analyzing all the ratios involved in the three 
BPT diagrams, we concluded that the evolution of galaxies in 
the three BPT diagrams is shown through the [O m] /!5007/Hy8 
ratio. Since this ratio is a metallicity indicator, any decrement 
in metallicity will result in higher values of the [O m] /15007/Hy6 
ratio. 

5. Evolution of the mass-metallicity and 
luminosity-metallicity relations 

It has been demonstrated that the metallicity and mass of SF 
galaxies are strongly correlated, with massive and luminous 
galaxies showing higher metallicities than less massive galaxies 
(see Sect. 1). 

The masses of our galaxies were estimated using the 
STARLIGHT code, which fit an observed spectrum with a com- 
bination of 150 SSPs from the evolutionary synthesis models of 
Bruzual & Chariot (2003), computed using a Chabrier (2003) 
initial mass function between 0.1 and 100 M , and "Padova 
1994" evolutionary tracks. The 150 base elements span 25 ages 
between 1 Myr and 18 Gyr, and six metallicities fromZ = 0.005 
to 2.5 Z . As argued by Mateus et al. (2006), the inclusion of 
very low Z SSPs in the base inevitably leads to larger stellar 
masses. A comparison with the Kauffmann et al. (2003b) mass 
estimates, which are based on a library of model galaxies con- 
structed with Z > 0.25 Z , results in systematic discrepancies of 
about 0.1 dex (for details see Mateus et al. 2006). The masses of 
our galaxies were corrected for aperture effects based on the dif- 
ferences between the total galaxy magnitude in the r band, and 
the magnitude inside the fiber, assuming that the mass-to-light 
ratio does not depend on the radius (see Mateus et al. 2006 for 
details). 

A histogram of our mass estimates is shown in Fig. 7, where 
a larger fraction of massive galaxies are observed at highest red- 
shifts. In Fig. 8 we derived the M - Z and L-Z relations for the 
galaxies of our sample. As explained in previous sections, galax- 
ies of the zo and z\ samples are complete in luminosity, while for 



the zi and Z3 samples, the completeness criterium is not taken 
into account. 

The M - Z relation of T04, which is valid over the range 
8.5 < log(M Jtor /M ) < 11.5, shows a steep M - Z relation for 
masses from 10 8 5 to 10 10 5 M that flattens at higher masses. In 
such study, T04 analyzed galaxies with redshift ranges 0.005 < 
z < 0.3. It is important to notice that in our M - Z relation, 
a flatness is not observed for masses > 10 10 5 for the redshift 
range zo, (see Fig. 8a), but this flatness is observed for the higher 
redshift samples. Then, the flatness observed by T04 depends on 
the redshift range observed. In order to establish the bias-free 
M - Z relation for local galaxies, Kewley & Ellison et al. (2008) 
recalibrated the M - Z relation of T04 with galaxies at 0.04 < z 
< 0.1, since 0.04 is the minimum redshift to avoid fibers effects 
(Kewley et al. 2005). 

In Fig. 8 (M - Z and L-Z relations) the metallicity decre- 
ment for the Z3 redshift sample discussed in our previous articles 
(Lara-Lopez et al. 2009a, b) is also evident. In order to fit our 
local M - Z relation, we estimated the mode of the metallic- 
ity of the galaxies in mass bins of 0.1 dex and fit them with 
a second order polynomial (y=ao + a\x + a2X 2 ), with «o = 
-0.467, a\ = 1.611,a2 = -0.067. We also fit a second or- 
der polynomial to our M - Z relation for galaxies at zi, with 
«o = -0.632, «i = 1.557, «2 = -0.063. All the fits are shown in 
Fig. 9. According to them, our M - Z relation for the galaxies 
at Z3 is ~ 0.2 dex lower compared to our local galaxy sample. 
Additionally, in Fig. 9 we compare the M - Z fits from litera- 
ture at different redshifts with our results. At z ~ 0.15 we rep- 
resent the calibration of T04, at z ~ 0.07 the T04 recalibration 
of Kewley & Ellison et al. (2008), as well as our fit for the local 
(zo) M — Z relation. At higher redshift, we represent our fit to 
the M - Z relation for galaxies at Z3, as well as the fit of Erb et 
al. (2006a) at z ~2.2 scaled to the T04 metallicity calibration. 
Due to their high redshift, Erb et al. (2006a) used the N2 method 
and the calibration of Pettini & Pagel (2004) to estimate their 
metallicities. We converted their N2 metallicities to the R23 cali- 
bration of T04 with the metallicity conversions given in Kewley 
& Ellison et al. (2008). Even with the dispersion of our local 
sample, our M - Z fit is a little lower, but in a good agreement 
with those of T04 and Kewley & Ellison et al. (2008). Since 
we are using the T04 calibration of the R23 method to estimate 
metallicities, the main differences with the fit of T04 are the red- 
shift ranges, as discussed above, and the mass estimates, since 
T04 and Kewley & Ellison et al. (2008) adopted a Kroupa et al. 
(2001) IMF, while we are using a Chabrier (2003) IMF. 

Although redshift ranges are different, the comparison of our 
fit for the z 3 sample with the Erb et al. (2006a) data at z~2.2, 
which also use a Chabrier IMF, are similar in 12+log(0/H) (see 
Fig. 9). As will be explained in the next sections, our Z3 sam- 
ple is conformed mainly by spiral galaxies, while the Erb et al. 
(2006a) sample corresponds to a mix of morphological types. 
A possible explanation for the high metallicities of Erb et al. 
(2006a), or to the lower metallicities of our sample, is given by 
Calura et al. (2009) who, using models that distinguish among 
different morphological types through the use of different infall, 
outflow and star formation, reproduce the mass-metallicity rela- 
tion in galaxies of all morphological types, taking as a reference 
the observational M - Z relations of Kewley & Ellison (2008), 
Savaglio et al. (2005), Erb et al. (2006a), and Maiolino et al. 
(2008). In his work Calura et al. (2009) predicts that at any red- 
shift, elliptical galaxies will present the highest stellar masses 
and the highest metallicities, whereas the irregulars are the least 
massive and metallic galaxies, being spiral galaxies at an inter- 
mediate stage. 
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This means that, being our Z3 sample composed only by spi- 
ral galaxies, our metallicities will be lower than if our sample 
were composed by a mix of morphological types. The observed 
metallicities and SFRs for the Erb et al. (2006a) sample at z~2.2, 
according to the study of Calura et al. (2009), indicate that their 
galaxies are likely to represent a morphological mix, partly com- 
posed of spirals (or proto-spirals) and partly of ellipticals (or 
proto-ellipticals). Calura et al. (2009) predicts for the Erb et al. 
(2006a) galaxies at z~2.2 lower metallicities (~0.3 dex) if the 
sample were composed by spiral galaxies. Then, the similarities 
in 12+log(0/H) of our M - Z relation at zs with the one of Erb 
et al. (2006a) at z~2.2, could be explained by the morphological 
selection in each case: lower metallicities in our M-Z relation 
can be addressed to the prominence of spiral galaxies, whereas 
larger metallicities in the M - Z relation of the Erb et al. (2006a) 
is a consequence of a sample formed by a mix of morphological 
types. In other words, the morphology of the galaxies is crucial 
in deriving and comparing the metallicity and the M-Z relation. 

An additional point to take into account, is that our samples 
plotted in Fig 9 are selected with differente magnitude complete- 
ness, then, our Z3 redshift sample is high luminous and massive 
than our local one. This must be taken into account when com- 
paring both M-Z relations. Unfortunately, our zq sample does 
not have a significant number of galaxies in the same absolute 
magnitude range of the Z3 sample, making imposible to generate 
a local M-Z relation comparable in luminosity to the Z3 one. 

Comparing our M-Z relation with that of high-z samples in 
the literature, we have the M-Z relation of Savaglio et al. (2005) 
at z ~ 0.7. In their study, they found only a slightly decrement in 
metallicity for galaxies at z ~ 0.7 compared with the local one of 
T04, which is inconsistent with the ~ 0.2 dex decrement found in 
our Z3 M-Z relation. Nevetheless, as pointed out by Rodrigues et 
al. (2008), the sample of Savaglio et al. (2005) have spectra with 
low S IN and spectral resolution, as well as extinction problems. 
With a more consistent result, Rodrigues et al. (2008) generated 
the M-Z relation for galaxies at z ~ 0.7, finding a decrement 
in metallicity of ~0.3 dex compared with the local one of T04. 
Unfortunately, given the small range in their stellar masses, it 
was not possible to constrain the evolution of the shape of their 
M-Z relation. 

There are two main ways to explain the origin of the M-Z 
relation. The first one is related to the well-known effect of 
downsizing (e.g. Cowie et al. 1996, Gavazzi & Scodeggio 1996), 
in which lower mass galaxies form their stars later and on 
longer time-scales than more massive systems, implying low 
star formation efficiencies in low-mass galaxies (Efstathiou 
2000; Brooks et al. 2007; Mouhcine et al. 2008; Tassis et al. 
2008; Scannapieco et al. 2008, Ellison et al. 2008). Therefore, 
low-mass galaxies are expected to show lower metallicities. 
Supporting this scenario, Calura et al. (2009) reproduced the 
M-Z relation with chemical evolution models for ellipticals, 
spirals and irregular galaxies, by means of an increasing effi- 
ciency of star formation with mass in galaxies of all morpho- 
logical types, without the need for outflows favoring the loss of 
metals in the less massive galaxies. In a recent study that sup- 
ports this result for massive galaxies, Vale Asari et al. (2009), 
model the time evolution of stellar metallicity using a closed- 
box chemical evolution picture. They suggest that the M-Z 
relation for galaxies in the mass range from 10 9 8 to 10 1L65 M 
is mainly driven by the star formation history and not by inflows 
or outflows. 

A second scenario to explain the M-Z relation is attributed 
to metal and baryon loss due to gas outflow, where low-mass 
galaxies eject large amounts of metal-enriched gas by super- 




Fig. 7. Normalized mass histogram for all our samples. Dark 
solid line represent galaxies at zo, dashed line galaxies at z\, 
point dashed line galaxies at Z3, and clear solid line galaxies at 
Z3- The histograms were constructed do not taking into account 
the completeness for galaxies at zi and Z3- 

novae winds before high metallicities are reached, while mas- 
sive galaxies have deeper gravitational potentials which helps to 
retain their gas, thus reaching higher metallicities (Larson 1974; 
Dekel & Silk 1986; MacLow & Ferrara 1999; Maier et. al. 2004; 
T04; De Lucia et al. 2004; Kobayashi et al. 2007; Finlantor & 
Dave 2008). As pointed out in the high-resolution simulations 
of Brooks et al. (2007), supernovae feedback plays a crucial role 
in lowering the star formation efficiency in low-mass galaxies. 
Without energy injection from supernovae to regulate the star 
formation, gas that remains in galaxies rapidly cools, forms stars, 
and increases its metallicity too early, producing a M-Z relation 
too flat compared to observations. 

An additional interpretation of the M-Z relation is linked 
to some properties of star formation, as the IMF. Koppen et al. 
(2007) suggested that the M-Z relation can be explained by a 
higher upper-mass cutoff in the IMF in more massive galaxies. 

Finally, we also generated the L-Z relation for our redshift 
samples (see Fig. 8b). Nevertheless, our Z3 sample is restricted 
to a small range in luminosity, making impossible to fit a curve. 
The local L-Z relation is well established by e.g. T04, then, 
due to our small luminosity range at Z3 we can not conclude any- 
thing about the L-Z relation. Nevertheless, the M-Z relation 
has demonstrated to be more stronger and tighter than the L-Z 
relation, confirming that stellar mass is a more meaningful phys- 
ical parameter than luminosity when both are compared with gas 
metallicity (Savaglio et al. 2005). 

5.1. Evolution of the mass-to-light ratio 

As explained by Erb et al. (2006a), at higher redshifts, the M-Z 
relation is clearly more physically meaningful than the L-Z 
relation. A corollary is that the local L-Z relation is simply a 
result of the strong correlation between mass and luminosity at 
low redshift. 

We also analyzed the evolution of the mass-to-light (M/L) 
ratio (see Fig. 10). For a given mass, we observed higher lumi- 
nosities for the Z3 sample compared with the local one, which 
means lower M/L ratios as redshift increase. In order to observe 
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Fig. 8. (a) Relation between the stellar mass and 12+log(0/H) (M - Z relation), (b) Relation between the absolute Petrosian r 
magnitude and 12+log(0/H) (L - Z relation) for our sample of galaxies. The cut observed in Fig. (b) for the zi sample, is due to 
the 5577 A sky line (see the text). In both relations, white contours represent, from outside to inside, 15, 30, 50, 70, and 90% of the 
maximun density value of the zo redshift sample (black dots), and are plotted only as a visual aid. Colors and symbols follow the 
same code used in Fig. 5. 




Fig. 9. Evolution of the mass-metallicity relation observed at 
different redshifts. Point-dashed line represents the curve of 
Tremonti et al. (2004) at z~0.15, dashed line represents the T04 
recalibration of Kewley & Ellison (2008) at z~0.07, circles are 
the mode metallicity in log(M s , flJ /M ) bins of 0. 1 for our zo sam- 
ple, solid dark and clear curve represents our fit for zo (mode 
bins) and Z3, respectively. Dotted line represents the fit of Erb et 
al. (2006a) at z~2.2. 



this evolution, we fit a line (y = ciq + a\x) to the zq and Z3 redshift 
samples (see Fig. 10), obtaining «o = -0.924 and a\ = -0.544 
for zo, and ao = -18.686 and a\ = -1.304 for Z3. The variation 
in M/L at a given rest frame optical luminosity can be as much 
as a factor of ~70 (Shapley et al. 2005), which means that for 
any range in luminosity there exist an extended range of stellar 



Fig. 10. Mass versus Petrosian absolute k-corrected magnitude 
for all our redshift samples. White contours represent, from out- 
side to inside, 10, 25, 50, and 80% of the maximun density value 
of the zo sample. Dashed line represents the fit to the zo sam- 
ple, while point-dashed line the fit to the Z3 sample. Colors and 
symbols follow the same code used in Fig. 5. 



masses. This large variation in M/L explains the lack of corre- 
lation in the L - Z relation for the Z3 sample compared to the 
local relation. For a small range of absolute magnitudes in the 
Z3 sample, we have a widely range of masses, making possible 
to generate a M - Z relation. At higher redshifts the effect is 
the same, as pointed out by Erb et al. (2006a), finding for star- 
forming galaxies at z~2.2 in a small range of luminosity, a wide 
range of stellar mass. 
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6. Morphology indicators and SFR 

The variation of SFR activity and young stellar content along 
the Hubble sequence is one of the most recognizable features of 
galaxies. In fact, this variation in stellar content is part of the ba- 
sis of the Hubble classification itself (Hubble 1926), and under- 
standing its physical nature and origin is fundamental to com- 
prehend the galaxy evolution (for a review see Kennicutt 1998 
and Kennicutt et al. 1994). The general picture, presented by 
Roberts (1963), Searle et al. (1973), Larson & Tinsley (1978), 
and Kennicutt et al. (1994), point that early-type galaxies (types 
SO-Sb) represent systems which formed most of their gas into 
stars on timescales much less than the Hubble time, while the 
disks of late-type systems (Sc-Im) have formed stars at roughly 
a constant rate since they formed. 

For most of the 20th century, catalogues of morphologi- 
cally classified galaxies were compiled by individuals or small 
teams of astronomers (e.g. Sandage 1961, de Vaucouleurs 1991). 
Nowadays selection criteria is based on galaxy properties such 
as color, concentration index, spectral features, surface bright- 
ness profile, other structural parameters, or some combination of 
these (e.g. Strateva et al. 2001, Abraham et al. 2003, Kauffmann 
et al. 2004, Conselice 2006, Scarlata et al. 2007). With the advent 
of modern surveys, such as the SDSS and with the participation 
of thousands of volunteers, it was possible the development of 
the Galaxy Zoo project (e.g. Lintott et al. 2008), providing vi- 
sual morphological classification for more than 10 7 galaxies. 

In this section, we will focus on the relations at different red- 
shifts between the SFR, metallicity, mass, and morphology of 
SF galaxies selected with the Kauffmann et al. (2001) criteria, as 
explained in the Section 2. 

6.1. Evolution of the SFR 

We estimate the SFR with the Ha emission line flux following 
the Kennicutt (1998) expression: 

SFR[M yr 1 ] = 7.9 x 10~ 42 L(Ha) [ergs s _1 ], (4) 

where L(Ha) denotes the intrinsic Ha luminosity, and Ha is cor- 
rected by dust extinction and underlying stellar absorption as ex- 
plained in Sec. 2. This calibration is derived from evolutionary 
synthesis models that assume solar metallicity and no dust, and 
is valid for a T c =10 4 K and case B recombination 

Recent studies have explored the relationship between the 
stellar mass and the SFR in galaxies at different redshifts. It has 
been shown that SFR critically depends on the galaxy mass both 
at low and high redshifts (e.g. Gavazzi et al. 2002, Brinchmann 
et al. 2004, Dickinson et al. 2004, Feulner et al. 2005, Papovich 
et al. 2006). 

In our sample, galaxies with high SFRs are more abundant 
at higher redshifts (see Fig. 1 1), a fact already observed in non 
biased samples (e. g. Noeske et al. 2007a). In Fig. 1 la, we show 
12+log(0/H) against log(SFR). Note that, although our Z3 sam- 
ple of galaxies is biased to the most luminous and massive galax- 
ies, the observed decrement of ~ 0.1 dex in 12+log(0/H) found 
in Lara-Lopez et al. (2009a, b) is also present. Regarding the zo 
sample of galaxies, there is a clear sequence with galaxies go- 
ing toward higher values of SFR as metallicity increases. This 
tendency can be explained from the zo sample in the M - Z rela- 
tion of Fig. 8, where massive galaxies corresponds to the highest 
metallicity galaxies, and for more massive galaxies, we expect 
higher SFRs (see Fig. 11). Also, we can slightly appreciate a 
population of galaxies with higher SFR (see Fig. 11a). This pop- 
ulation will form a tail when the mass is taken into account, 



as will be shown in Fig. 1 lb. As redshift increases, we appre- 
ciate in Fig. 11a for galaxies at zi and Z3, a flattening of the 
SFR vs. 12+log(0/H) relation, with most of the galaxies show- 
ing log(SFR) between 1 and 2. 

In Fig. lib, we show the log(M Jtor /M ) versus log(SFR) 
plot. Galaxies at zo show a main sequence, where massive galax- 
ies have higher SFRs. This main sequence was identified by 
Noeske et al. (2007a), studying galaxies with redshifts from 
0.2 to 1.1, finding that this main sequence moves as a whole 
to higher SFR as redshift increases. The SSFR, defined as the 
total SFR divided by the stellar mass, reflects the strength of 
the current burst of star formation relative to the underlying 
galaxy mass. Deep galaxy surveys have consistently found that 
the SSFR depends strongly on both M and redshift, with the 
bulk of star formation occurring earlier in massive galaxies than 
in less massive systems (Guzman et al. 1997, Brinchmann & 
Ellis 2000, Juneau et al. 2005, Bauer et al. 2005, Bell et al. 2005, 
Perez-Gonzalez et al. 2005, Feulner et al. 2005, Papovich et al. 
2006, Caputi et al. 2006, Reddy et al. 2006). 

We also analyzed the evolution of the SSFR as a func- 
tion of the metallicity and the stellar mass. Interestingly, in the 
12+log(0/H) vs. SSFR diagram (Fig. 12a), the observed popu- 
lation of Fig. 1 la at zo is more evident, showing a higher SSFR 
(SSFR > -10) than the other galaxies at the same redshift. We 
are going to investigate this tail in more detail in the next sub- 
section. 

In Fig. 12b, we show the log(M SMr /M ) versus log(SSFR) 
plot. The SSFR increase with redshift, showing for more mas- 
sive galaxies a tendency toward lower SSFR values, which is in 
agreement with the results of Noeske et al. (2007b), for galaxies 
with z > 0.2. Massive galaxies shows lower SSFR because they 
probably have low gas fractions and have thus nearly finished as- 
sembling their stellar mass (Erb et al. 2006b, Reddy et al. 2006). 
On the other hand, the presence of dust has demonstrated to play 
an important role deriving the SSFR, as shows by Pannella et 
al. (2009), getting for dust free galaxies a flat SSFR for redshift 
bins centered at z ~ 1.6 and 2.1, instead of a drop with increasing 
mass. In our samples, the SFR and SSFR were derived through 
the dust corrected Ha flux, showing that for redshifts less than 
0.4, this flatness is absent. 

We found also evidences of two populations in Fig. 12b, one 
at redshift zo which concentrate in a square delimited by the con- 
tour plots, and another one showing SSFR > - 10, which is more 
evident for the higher redshifts samples, but it is also evident in 
the zo sample. This population of galaxies will be explained in 
the next subsection as a consequence of the morphology of the 
galaxies. 

6.2. The SSFR as a morphology indicator 

In this section, we investigate the morphology of the galaxies in 
our sample with the aim of clarify if the tail observed in Figs. 
1 1 and 12, with higher SFR and SSFR in the zo sample, respec- 
tively, is related to specific morphological types. We will focus 
on samples zo and zj, because those samples show their galax- 
ies uniformly distributed on mass. Galaxies of samples z\ and zi 
show systematic problems due to sky lines, and to incomplete- 
ness, respectively. 

We used the g - r color, the concentration index c=R9o/Rso 
(e.g. Park &Choi et al. 2005), and the SSFR (e.g. Salim et al. 
2009), which are the most common indexes to segregate early 
from late-type galaxies (see Fig. 13). We used k-corrected fiber 
colors for all the galaxies samples. 




Fig. 11. Metallicity and mass versus log(SFR), contours correspond to the zo sample in both plots. White contours represent, from 
outside to inside in panel a: 15, 30, 50, 70, and 90%, and in panel b: 5, 15, 35, 65, and 85% of the maximun density value of the zo 
sample. Colors and symbols follow the same code used in Fig. 5. 
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Fig. 12. Metallicity and mass versus log(SSFR), contours correspond to the zo sample in both plots. White contours represent, from 
outside to inside in panel a: 10, 20, 40, 60, and 90%, and in panel b: 5, 15, 35, 60, and 80% of the maximun density value of the zo 
sample. Colors and symbols follow the same code used in Fig. 5. 



The color index has been commonly used as an early and 
late type morphological classificator (Baldry et al. 2004, Faber 
et al. 2007, Wang et al. 2007, Lee et al. 2007). Strateva et al. 
(2001) found that the integrated observed frame u-r shows a bi- 
modal distribution, however, they have shown that, when divided 
at u — r = 2.22, the early and late type subsets have significant 
contamination, reaching about 30% for a sample with visually 
identified morphological types. Because the u band shows large 
errors for the SDSS galaxies, we decided to use g - r < 0.6 
(e.g. Schawinski et al. 2009, Masters et al. 2009), which allows 
separate early from late-type galaxies. As observed in Fig. 13a, 



~97% of the Z3 sample of galaxies correspond to a late-type 
morphology. 

The concentration index c=R9o/R5o has been successfully 
used in segregating late (c < 2.86) from early-type (c > 2.86) 
subsets (e.g. Shimasaku et al. 2001, Strateva et al. 2001, Goto et 
al. 2002, Nakamura et al. 2003, Deng et al. 2009). Nevertheless, 
contamination in the early and late-type subsets separated using 
the concentration index, is typically about 20% (Yamauchi et al. 
2005, Shimasaku et al. 2001). Using both, color and concentra- 
tion index, Park & Choi et al. (2005) used the color-color space 
u-r versus A(g - r) and the concentration index c~'=R5o/R9o 
~ 0.35 as a reliable morphological classificator. In Fig. 13b, we 




Fig. 13. Principal morphological classificators. In panel a, b, and c, we show the log(M s , af /M ), concentration index c=R9o/Rso, and 
log(SSFR) versus g — r color, respectively. Black dots and yellow circles represent galaxies at zo and Z3, respectively. In each panel 
the dashed line shows the standard limit to segregate early from late-type galaxies, with the zj, galaxies concentrated in the late-type 
region in each panel. 
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Fig. 14. Metallicity and Mass versus log(SSFR). Black dots and yellow circles represent galaxies at zo and Z3, respectively, while 
green dots, represent late-type galaxies at zo selected with log(SSFR) > -10. 



show the concentration index cvs.^-r for galaxies at zo and Z3, 
with the ~ 95% of the Z3 galaxies corresponding to late-type (c 
< 2.86) galaxies. 

Finally, the SSFR has been used as an indicator of early 
and late-type morphology (e.g. Wolf et al. 2009, Salim et al. 
2009) since late-types have blue colors and high SSFRs, while 
early-types have red colors and low SSFRs. For our Z3 sample, 
~ 89% of the galaxies have log(SSFR) > -10 (see Fig. 13c), 
and as reported by Salim et al. (2009), blue actively star-forming 
galaxies has log(SSFR) > -10, while lower values would corre- 
spond to the green valley and red-sequence galaxies. As argued 
by Weinmann et al. (2006), the use of the SSFR would give us 
important clues in determining the morphological galaxy type, 
since for example, a genuine SF disk galaxy may appear red 
due to strong extinction (e.g. when seen edge-on), and thus be 
classified as early-type based on its color, while the SFR and 
morphology quantifiers would classify it as a late-type galaxy. 

Our sample of galaxies at Z3 is mainly composed by late-type 
galaxies, as indicated by any of the morphological classificators 
discussed above. We decided to take these galaxies as a refer- 
ence to delimitate the late-type zone. After trying all the dis- 



cussed methods, in order to obtain a reliable morphological clas- 
sification, and following the location of the Z3 galaxies in Fig. 
13c, we conclude that the best method is to use both, log(SSFR) 
> -10, and a color g - r < 0.6. After applying this criterion to 
our zo sample, we end with 7967 galaxies classified as late-type, 
corresponding to the ~13% of the original sample. If we repeat 
the comparative analysis between SFR and metallicity, and SFR 
against stellar mass (see Fig. 14) for galaxies at zo and z\, the 
tail with higher SFR identified in previous sections corresponds 
to late-type galaxies. Note that the separation between late and 
early-type in the stellar mass versus SFR is a straight frontier 
(see Fig. 14b), because the separation criterion is the SSFR. 

In order to alternatively assess this statement, we com- 
pare our results with selected mock galaxy samples from the 
Millennium simulations (Springel et al. 2005). We used the 
Bower2006a catalog (Bower et al. 2006), which give us red- 
shift, SDSS k-corrected colors, stellar mass, and Ha luminosity, 
among other information. From the original catalog, we selected 
randomly 35000 galaxies at z=0 and z=0.4 with the Ha line in 
emission in order to be able to compare with our galaxies. As 
shown in Fig. 15, there is a clean separation between a red se- 
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Fig. 15. Mass versus g — r color for Millennium galaxies. Black 
and yellow dots represent galaxies at zo and zj, respectively. 

quence, formed by galaxies with g—r > 0.8, the green valley with 
0.6 < g - r < 0.8, and the blue cloud with g - r < 0.6. It can be 
observed how galaxies move towards late-type galaxies as red- 
shift increases. Then, this results are consistent with observing 
mainly late-type galaxies at redshift Z3- 

7. Summary and Conclusions 

We analyzed a sample of emission line galaxies selected in four 
redshift intervals from ~0 to 0.4 in bins of 0.1, taking into ac- 
count the magnitude completeness of every redshift interval. In 
this paper we introduced the S2N2 diagram as a star-forming, 
composite, and AGNs galaxy classificator, we estimated metal- 
licities using the R23 method and analyzed the evolutive effects 
of galaxies from the three BPT diagrams. Additionally, we stud- 
ied the evolution of the M - Z and L-Z relations, and ana- 
lyzed the evolution and implications of the galaxy morphology 
in the SFR-mass and metallicity relations. From these analysis 
we conclude the following: 

- Using the KewOl photoionization grids, and the Kauf03 and 
KewOl SF, and starburst limit respectively, in the [Nn] /Ha 
vs [O m] /15007/Hy6 diagram, we have demonstrated that the 
S2N2 is a well-behaved diagnostic diagram efficiently clas- 
sifying star-forming, composite, and AGNs galaxies. 

- We analyzed the galaxy evolution using the three main BPT 
diagrams: [N 11] /Ha, [S 11] /Ha, and [0 1] ^6300/Ha vs [O m] 
/15007/H/? in our four redshift bins, observing an evolution 
toward higher values of the [O m] A5007fH/3 ratio. This evo- 
lution is a consequence of the metallicity evolution as red- 
shift increases, reflected in the three BPT diagrams, because 
the ratio [O m] /15007/H/3 is a good metallicity indicator. As 
a result, a metallicity decrement will be reflected in higher 
values of this ratio. 

- We analyzed the evolution of the M - Z and L-Z relations, 
observing that at higher redshift values, both relations evolve 
towards lower values of metallicity. We discovered that the 
flat zone of the M - Z relation reported by Tremonti et al. 
(2004) for galaxies with log(M J/flr /M ) > 10.5, is mainly 
constituted by galaxies at z > 0.1 (samples at zi, Zi and Z3). 
Galaxies at zo redshift could be fitted with a linear function. 
Our M - Z relation at redshift Z3 is ~0.2 dex lower than our 
local one. 



- Our fit to the M - Z relation for sample Z3 is in agreement 
with the one of Erb et al, (2006a) at z ~ 2.2. We attribute 
this similarity to the galaxy morphology in the different sam- 
ples, since our Z3 sample is conformed by late-type galaxies, 
while the sample of Erb et al. is composed by a mix of early 
and late-type galaxies. According to Calura et al. (2009), the 
M — Z relation of late-type galaxies will have systematically 
lower metallicities than a M - Z relation conformed by a mix 
of early and late-type galaxies. 

- We analyzed the evolution of the mass-to-light ratio, observ- 
ing lower M/L ratios as redshift increase. For a small range 
of absolute magnitudes in the Z3 sample, we have a wide 
range of mass, making it possible to generate the M - Z re- 
lation, but difficult to generate the L-Z relation. 

- The decrement in metallicity observed in previous papers for 
galaxies at redshift Z3 (Lara-Lopez et al. 2009a,b) is also ob- 
served, even though in this study we are not restricting our 
galaxy luminosities as in our previous studies. 

- We estimated the SFR and SSFR for our sample of 
galaxies and analyzed its relation with 12+log(0/H) and 
log(M J/£lr /M ), confirming the existence of a main sequence 
reported by Noeske et al. (2007) in the log(SFR) vs. 
log(M J , fl) /M ) plot. Consistently, we found that higher SFRs 
and SSFRs increase with redshift. 

- We analyzed the morphology of our galaxies through the 
g — r color, the concentration index R90/R50, and the SSFR, 
concluding that the best method to determine the morphol- 
ogy was combining both, a color of g - r < 0.6, and a 
log(SSFR)>10 for selecting late-type galaxies. 

- Our Z3 sample of galaxies is mainly formed by late-type 
galaxies, a fact that helped us to classify morphological types 
at lower redshift. The fact that at higher redshift the fraction 
of late-type galaxies is larger, was confirmed by using mock 
galaxy catalogues from Millennium simulations. 

- We found at the higher redshift, a population with higher 
SFR and SSFR than the galaxies in the zo sample. After clas- 
sifying late and early-type galaxies in the zo sample, we real- 
ized that the observed tail showing higher SFR and SSFR is 
formed by late-type galaxies, demonstrating the connection 
of the galaxy morphology with the SFR in a new fashion. 

Our work provide a useful tool for classifying galaxies with 
the S2N2 diagram, and demonstrating how galaxies evolve on 
the BPT diagrams as a consequence of metallicity evolution. 
We also analyzed the mass, metallicity and SFR relations, not- 
ing that galaxies in the redshift sample Z3 have lower values of 
metallicity, higher SFRs, and morphology indicators associated 
to late-types. In this study we pointed out the importance of the 
morphology of galaxies when deriving conclusions. Since a sam- 
ple conformed by late-type galaxies will show lower values of 
metallicity than ones formed by a mix of morphological types. 
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